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Abstract 

Background: Insiglits into tlie meclianisms associated witli cliemotlierapy-resistance are important for 
implementation of tlierapeutic strategies and for unraveling the mode of action of chemotherapeutics. Although 
cidofovir (CDV) has proven efficacious in the treatment of human papillomavirus (HPV)-induced proliferation, no 
studies concerning the development of resistance to CDV in HPV-positive tumor cells have been performed yet. 

Methods: From the cervical carcinoma SiHa cells {SiHdparenta^ which are HPV-16 positive, cidofovir-resistant cells 
(SiHacDi/) were selected, and differential gene expression profiles were analyzed by means of microarrays. We 
examined in vitro phenotyping of resistant cells compared to parental cells as well as tumorigenicity and 
pathogenicity in a mouse-xenograft model. 

Results: SiHacoi/ had a resistant phenotype and a reduced growth both in vitro and in vivo. A markedly diminished 
inflammatory response (as measured by production of host- and tumor-derived cytokines and number of 
neutrophils and macrophages in spleen) was induced by SiHacoi/ than by S\[^3parentai in the xenograft model. 
Gene expression profiling identified several genes with differential expression upon acquisition of CDV-resistance 
and pointed to a diminished induction of inflammatory response in SiHacoi/ compared to S\h3parentai- 

Conclusions: Our results indicate that acquisition of resistance to cidofovir in SiHa cells is linked to reduced 
pathogenicity. The present study contributes to our understanding on the antiproliferative effects of CDV and on 
the mechanisms involved, the inflammatory response playing a central role. 
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Background 

Three acyclic nucleoside phosphonate analogues (ANPs), 
i.e. tenofovir (PMPA), adefovir (PMEA) and cidofovir 
(CDV), are approved for the treatment of viral infections 
[1,2]. Tenofovir and adefovir are active against retrovi- 
ruses and hepadnaviruses, their oral prodrug forms 
being approved for therapy of HIV (PMPA) and of 
chronic hepatitis B virus infections (PMPA and PMEA). 
Although CDV is formally licensed for treatment of 
cytomegalovirus retinitis in AIDS patients, it is often 
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used off-label for the management of diseases caused by 
several DNA viruses, including adeno-, pox-, papilloma-, 
polyoma-, and herpesviruses others than cytomegalo- 
virus [3-6]. 

Besides their well-recognized antiviral properties, some 
ANPs have shown anticancer potency. For instance, 
PMEA, PMEDAP, PMEG, and prodrugs of PMEG [i.e. 
cPr-PMDEDAP, GS-9219 and GS-9191] showed marked 
cytotoxic properties in vitro [7-9]. Additionally, in vivo 
antitumor activities for these compounds have been 
described in different animal models: GS-9219 in a pet 
dog model of non-Hodgkin's lymphoma [10] and cPr- 
PMEDAP in a rat choriocarcinoma tumor model [11]. A 
close correlation between the cytostatic activities of 
PME derivatives and the inhibitory effects of their active 
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metabolites (diphosphate forms) on cellular DNA poly- 
merases a, 5, and e has been established. In these stud- 
ies, PMEG- diphosphate (PMEGpp) emerged as the most 
potent chain-terminating inhibitor of cellular DNA poly- 
merases [12,13]. The utility of PMEG as an anticancer 
agent is limited by poor cellular permeability and tox- 
icity [13,14] and prodrugs, such as GS-9191 and GS- 
9219, were designed to increase the permeability and 
accumulation of PMEGpp in the cells [10,13]. 

Cidofovir represents also an ANP with marked anti- 
proliferative effects but unlike PMEG, the effects of 
CDV-diphosphate (CDVpp) on cellular DNA polymeri- 
zation are weak [inhibition constant (Ki) of CDVpp for 
cellular DNA polymerase-a of 51 [iM versus 0.55 [iM for 
PMEGpp]. In addition, CDVpp is not an obligate chain 
terminator [12,13] and, in contrast to PMEG, CDV has 
been used to manage human papillomavirus (HP V) -in- 
duced benign and malignant hyperproliferation with 
minimal if any side-effects, as described in several case 
reports and some phase II/III clinical trials [15-20]. Re- 
cently, a phase II clinical trial was conducted in Belgium 
to evaluate the safety and efficacy of CDV in the treat- 
ment of high grade cervical lesions (NCT01303328). Full 
data analysis of this Phase II clinical trial will be pro- 
vided during the next months. 

Cidofovir antitumor properties were also demonstrated 
in different animal models of tumors related to trans- 
forming viruses, including Epstein-Barr virus-associated 
nasopharyngeal carcinoma [21] and HP V- induced cer- 
vical carcinoma [22-24] xenografts in athymic-nude 
mice, polyomavirus-induced hemangiomas in rats [25] 
and hemangiosarcoma development in mice [26]. Also, 
CDV proved effective against cottontail rabbit papillo- 
mavirus in the domestic rabbit model [27]. 

We have recently shown that besides inhibition of 
tumor growth, intratumoral CDV administration had a 
beneficial effect on the pathology associated with the 
growth of cervical carcinoma cells in athymic nude mice 
as demonstrated by a favorable effect on body weight gain, 
reduced splenomegaly and lower inflammatory state in an- 
imals that received the compound versus the placebo- 
treated group [24]. Furthermore, a whole genome gene 
expression profiling performed on CDV-treated malignant 
cells and normal keratinocytes allowed us to identify 
unique signatures in tumor cells compared to normal ker- 
atinocytes pointing to a selective drug effect [28]. Among 
the functions that were distinctly regulated by CDV in ma- 
lignant and normal cells, the acute phase response was 
found exclusively activated in transformed cells but not in 
normal keratinocytes. In addition, cell cycle regulation and 
DNA repair by homologous recombination was only acti- 
vated in normal cells [28] . 

There are several mechanisms by which cancer cells 
develop drug-resistance and this is often a multi- 



factorial process. Understanding the mechanisms leading 
to development of drug-resistance is crucial for the 
implementation of therapeutic strategies, for providing 
insights into the effects of anticancer drugs on specific 
cellular functions, and also for predicting how acquisi- 
tion of drug-resistance impacts tumor igenicity and pa- 
thogenicity. Therefore, we established, from the cervical 
carcinoma cell line SiHa (HPV16^), a CDV-resistant cell 
subline (denoted SiH^lcdv) by stepwise dose escalation 
of CDV. We investigated the in vitro and in vivo 
phenotyping and growth rate of SiHac^y compared to 
parental cells {SiHsLparentai)- Also, we evaluated the differ- 
ential gene expression profiles between SiHsiparentai and 
SiHac^y by microarray analysis in order to identify 
genes changing expression upon selection of cells for 
CDV-resistance. In the present study, we focused on the 
analysis of functions and pathways involved in the in- 
flammatory response that changed in SiHa cells follow- 
ing acquisition of CDV-resistance. Importantly, we also 
examined whether SiHa cells that acquired resistance to 
CDV were impaired in pathogenicity in the xenograft 
model. 

Results 

In vitro phenotyping of SiHacoi/ 

SiHa cells {SiHsiparentai) were selected for CDV-resistance 
{SiHdiCDv) following continuous in vitro exposure to 
the drug for approximately 45 passages. The resulting 
SiHacDK presented a reduced growth rate compared to 
SiHdiparentai (doubling time of 40 h versus 26 h) (Figure 1). 

A stable CDV-resistant phenotype in the absence 
of selective drug pressure was found for the SiH^lcdV' 
When evaluated in terms of cell growth inhibition 




Days 

Figure 1 Growth rate of SiHacD\/ compared to SiHaparentai- Cells 
were seeded at a density of 2.5 x 10^ cells per well in a volume of 
0.1 ml into 96-well microtiter plates and allowed to proliferate. At 
several time-points, cells were trypsinized and the number of cells 
was determined with a Coulter Counter. Doubling time (DT) was 
calculated with the formula: DT= (t - to)/(log2N - log2No), where 
t and to are the times at which the cells were counted, and N and 
No are the cell numbers at times t and to. 

V J 
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(Table 1), a fold-resistance of >100 against CDV was 
determined after 7 days of incubation with the drug. 
SiHacDK displayed < 10-fold cross-resistance to the 
cytosine analogue Ara-C and to two unrelated ANPs 
(i.e. PMEG and cPr-PMEDAP). It should be noted that 
the antiproliferative effects of CDV (but not of PMEG 
or cPr-PMEDAP) for SiHsiparentai were time-dependent, 
pointing to a different mechanism of antiproliferative 
effects for these drugs, in agreement with our previous 
report [29]. 

An inhibition of 93% (SiHsiparentai) and 11% (SiHac^y) 
in the number of cells was afforded by CDV treatment 
at 158.7 [iM for 7 days. To compare CDV effects on in- 
duction of apoptosis in these cell cultures, annexin V 
and PI staining was performed. Annexin V stains phos- 
phatidylserine, a negatively charged phospholipid that is 
translocated from the inner leaflet of the plasma mem- 
brane to the outer leaflet during early apoptosis. Since 
PI does not enter into cells with intact membranes, it 
was used to identify necrotic cells. SiHsiparentai treated 
with CDV (at 158.7 [iM and 63.5 [iM) for 7 days showed 
increased percentage of apoptotic cells and diminished 
amounts of viable cells (Table 2). In contrast, SiHac^y 
were totally refractory to CDV-induced apoptosis while 
they were still able to respond to PMEG, albeit to a 
lower extent compared to SiH^Lparentai- No signs of cell- 
death by necrosis were seen in any of the two cell lines 
following treatment with either CDV or PMEG. 

Differentially expressed genes upon acquisition of CDV- 
resistance 

Gene expression profiling by microarray was performed 
to identify potential mechanisms associated with CDV- 
resistance. A total of 1,340 expression changes were 
identified in SiHslcdv compared with SiHsiparentah 777 
genes being upregulated and 563 downregulated. To val- 
idate the microarray results, transcript levels of four 



genes (randomly selected before bioinformatics analysis) 
were evaluated by qPCR, the expression patterns totally 
matching the microarray data (Additional file 1). 

Functional classification of differentially expressed 
genes showed that they were implicated in a variety of 
diverse and widely distributed functional categories and 
biochemical pathways: 12 functional categories and 106 
canonical pathways were associated with acquisition of 
CDV-resistance in SiHa cells (data not shown). 

In the present study, we focused on the inflammatory 
response (one of the 12 functional categories associated 
with acquisition of CDV-resistance), based on our previ- 
ous findings showing that CDV treatment of three ma- 
lignant cells (including SiHsiparentai) and primary human 
keratinocytes allowed the identification of acute phase 
response signaling' as a pathway exclusively modulated 
by CDV in transformed cells but not in normal cells 
[28]. When analyzing the immune response functional 
category in SiHac^y versus SiHsiparentah CDV-resistance 
was linked to a decrease (negative z-score) in four func- 
tional annotations: 'inflammatory response', activation of 
granulocytes', 'inflammation of organ', and activation of 
neutrophils' (Additional file 2). 

Twenty-one out of 106 canonical pathways affected 
by the changes in gene expression when comparing 
SiHdiCDv and SiHsiparentai Were related to immune re- 
sponse (Additional file 3). Several interleukin signaling 
pathways (such as lL-1 signaling', lL-6 signaling', lL-8 
signaling', lL-9 signaling', lL-10 signaling') as well as 
'interferon signaling', the endogenous danger signaling 
pathway 'HMGBl signaling', the prototypical proin- 
flammatory signaling pathway 'NF-kB signaling', the 
acute phase response signaUng', 'Toll-like receptor 
signaling', and 'MSP-RON signaling pathway', were 
among the several inflammatory response-related path- 
ways altered following acquisition of CDV-resistance in 
SiHa cells. 



Table 1 Growth inhibition of S'lHaparentai 


and S'lHacDv 


after incubation with CDV, PMEG, cPr-PMEDAP or Ara-C 


Compound 


Cell line 


CCso (MM) 






Fold-resistance 






Day 3 


Day 5 


Day 7 


at day 7 


CDV 


SiHdparental 


77.1 ±22.5.1 


15.2 ±13.6 


6.3 ± 6.0 


>100 




S\[-\3cDV 


>634.7 


>634.7 


>634.7 




PMEG 


SiHapQ^enra/ 


0.85 ± 1.14 


0.16±0.16 


0.10 ±0.07 


2.3 




SiHacDi^ 


1 .79 ± 2.34 


0.39 ± 0.33 


0.23 ±0.016 




cPr-PMEDAP 


S'\[^3parental 


1 .22 ± 0.67 


0.21 ±0.15 


0.21 ±0.06 


2.9 




SiHacDi/ 


3.44 ± 0.52 


1 .22 ± 0.34 


0.61 ±0.09 




Ara-C 


SiHdparental 


0.08 ± 0.04 


0.06 ± 0.02 


0.1 3 ±0.02 


8.8 




SiHacD\/ 


1.15 ±0.37 


0.99 ± 0.37 


1.15 ±0.49 





Cells were seeded at a density of 2.5 x 10 cells per well In a volume of 0.1 ml into 96-well microtiter plates and allowed to proliferate for 24 h. At this time, media 
containing different drug concentrations (tested in duplicate) were added (100 |al/well). After 3, 5, and 7 days of incubation, the cells were trypsinized and the 
number of cells was determined with a Coulter Counter. Antiproliferative effects were expressed as CC50 (50% cytostatic concentration), or concentration required 
to reduce cell growth by 50% (relative to the number of cells in the untreated control cell cultures). Data are mean ± standard deviation of at least two independent 
experiments. Fold-resistance was calculated as the ratio CC50 for SiHacD\/to CC50 for S\Haparentah 
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Table 2 Apoptosis detection by annexin V binding/PI uptake assay in SiHaparentai and SiHacov 



SlHaparental S'lHacDV 





Untreated 


CDV 158.7 mM 


CDV 63.5 mM 


PMEG 6.5 mM 


Untreated 


CDV 158.7 mM 


CDV 63.5 mM 


PMEG 6.5 mM 


Early apoptotic 


1 .2 ± 0.9 


22.1 ±0.6 


8.4 ± 1 .6 


64.2 ± 1 .2 


4.5 ± 2.0 


2.3 ±0.1 


2.4 ± 0.2 


30.6 ±4.3 


Late apoptotic 


0.6 ± 0.2 


5.0 ±5.6 


2.7 ±2.5 


6.0 ± 4.2 


0.9 ± 0.4 


0.6 ±0.1 


0.6 ±0.1 


3.9 ±0.7 


Viable 


97.2 ± 2.2 


72.2 ± 6.6 


88.0 ± 3.9 


29.4 ± 3.0 


94.5 ± 2.2 


96.3 ± 0.5 


96.5 ± 0.3 


64.9 ± 1 .2 


Necrotic 


0.4 ± 0.3 


0.7 ± 0.3 


0.9 ± 0.2 


0.4 ±0.1 


0.2 ±0.1 


0.8 ± 0.6 


0.5 ± 0.2 


0.7 ± 0.3 



Cells were seeded in 6-well microtlter plates at a density of 5 x 10"^ cells/well and allowed to proliferate for 24 h. Culture medium was removed and medium con- 
taining the compounds was added. At 7 days post-treatment, cells were harvested with EDTA/PBS (760 mg/L). After rinsing the cells with PBS, the cell pellets were 
simultaneously stained with annexin-FITC and PI using the annexin-V-FITC staining kit (BD Biosciences, San Jose CA) and analyzed by flow cytometry. Dual fluores- 
cence dot plots were obtained following combined annexin V-FITC and PI staining. Data represent the percentages in each quadrant of the dot plots and are the 
mean ± SD of at least two stainings. Viable (annexin VVPI") early apoptotic (annexin VVPI"), necrotic (annexin VVPr), and late apoptotic (annexin VVPI"^) cells. 



Using stringent criteria for microarray analysis (> 2- 
fold change in expression and p-value < 0.05), 173 genes 
related to the 'inflammatory response' function were 
identified as differentially expressed when comparing 
SiHsiCDv and SiHsiparentah 86 being upregulated and 87 
downregulated (Additional file 4). In order to visualize 
the interactions of differentially expressed genes involved 
in the inflammatory response, a network was con- 
structed based on the differentially expressed genes in- 
volved in the inflammatory response function (Figure 2). 
Gene networks represent intermolecular connections 
among interacting genes based on functional knowledge 
inputs stored in the IPKB (ingenuity Pathways Know- 
ledge Base). The different pathways and functions in 
the inflammatory response associated with CDV- 
resistance are indicated in this network. Among several 
genes changing expression upon acquisition of CDV- 
resistance, decreased gene expression of TGFBl (trans- 
forming growth factor beta-1), STATS [signal transducer 
and activator of transcription 3 (acute-phase response 
factor)], S0CS3 (suppressor of cytokine signaling 3), 
FOS (proto-oncogene c-Fos or activator protein 1), TLR3 
and TLR4 (toll-like receptor 3 and 4) and increased 
gene expression of CCNDl (cyclin D), CXCL2 [chemo- 
kine (C-X-C motif) ligand 1], CEBPB (interleukin 6- 
dependent DNA-binding protein) and STATl (signal 
transducer and activator of transcription- 1) appeared to 
play a central role in the changes in the inflammatory 
response that accompanied the development of CDV- 
resistance (Figure 2). These genes were implicated in 
several pathways and functions changed in SiHac^y versus 

In vivo phenotyping, tumorigenicity and pathogenicity of 

In a following step, we examined whether SiHac^y 
presented a resistant phenotype in vivo. Four weeks 
intratumoral treatment with CDV of athymic nude mice 
bearing the SiHac^y xenografts resulted in a moderate 
but not significant effect on tumor growth (Additional 
file 5). In contrast, the same treatment given to mice 



harboring the SiH3.parentai xenograft caused a remarkable 
and significant suppression of tumor growth, in agree- 
ment with our previous report [24]. Interestingly, tumor 
size was significantly lower in the SiHac^y cohort than 
in the SiHsiparentai group. 

Therefore, we investigated the kinetics of growth and 
the pathogenicity of SiH^lcdv in the xenograft model 
in athymic nude mice. The growth rate of SiH^lcdv tu- 
mors was significantly reduced compared to SiHsiparentai 
(Figure 3A). SiHdiparentai tumor size (351.6 ± 259.8 mm^) 
at week 3 was equivalent to that of SiH^lcdv (342.0 ± 
182.3 mm^) at week 5 (Figure 3A). In contrast to ani- 
mals inoculated with SiHsiparentah the cohort bearing 
SiHdiCDv tumors did not develop signs of wasting syn- 
drome. Body weight gain determined after subtraction of 
the tumor weight showed that the SiHsiparentai group 
gained significant less weight than tumor-free control 
animals (Figure 3B). Despite no statistically significant 
differences, mice with SiHac^y xenografts presented an 
intermediate body weight gain. 

Blood test analysis at week 5 showed similar results in 
the SiHdicDv tumor and tumor-free groups, except for a 
significant increase in the platelet count in mice bearing 
SiHac^y tumors (Table 3). By contrast, the SiHsiparentai 
group showed significant changes in several blood pa- 
rameters compared to either tumor-free or SiH^lcdv 
tumor mice. Hence, animals harboring SiHsiparentai tu- 
mors showed markedly higher white blood cell (WBC) 
count and lower RBC (red blood cell) count than tumor- 
free animals, most likely as a consequence of the growth 
of the xenograft and the subsequent induced inflamma- 
tion. The increase in the platelet number observed in 
the SiHdiparentai gfoup (but not in the SiHac^y cohort) 
was accompanied by an increase in the mean platelet 
volume. The growth of SiHsiparentai in mice also caused a 
marked decrease in hemoglobin concentration and in 
the mean corpuscular hemoglobin concentration, point- 
ing to anemia. Also, the liver enzymes aspartate amino- 
transferase and gamma-glutamyl transferase as well as 
lactate dehydrogenase were increased in the blood of the 
SiHdiparentai cohort. Elevation of creatine kinase (which is 



De Schutter et al. Molecular Cancer 2013, 12:158 
http://www.nnolecular-cancer.conn/content/1 2/1/1 58 



Page 5 of 16 



l\^YL'K 



CP: Clathrin-mediated Endocytosis Signaling 



B Cell Receptor Signaling 




CP: IL-9 Signaling 



CP: Interferon Signaling 



Figure 2 Network of inflammatory response and corresponding transcripts. Networks were constructed with IPA software using genes DE 
involved in 'inflammatory response'. A network is a graphical representation of the molecular relationships between molecules (nodes). The 
biological relationship between two nodes is represented as an edge (line) connecting two nodes. All edges are supported by information from 
the literature stored in the Ingenuity Pathways Knowledge Base. A solid line represents a direct interaction between two gene products and a 
dotted line means there is an indirect interaction. The intensity of the node color indicates the degree of up-regulation (red) or down-regulation 
(green) when comparing SiHacoi/ versus S\^3parentai cells. Canonical pathways identified by IPA in the networks are shown in blue while the 
functional annotations are shown in yellow. 



an indication of damage to muscle) and of creatinine 
(a sign of impaired Iddney function) was also seen in the 
SiWdiparentai gfoup. Taken together, these data clearly sig- 
nified a pronounced alteration of hematological and 
chemical blood parameters in the SiHdiparentai group and 
a lesser pathogenic effect of SiHac^y versus SiHsiparentai 
in the xenograft model. 



The SiWdiparentai cohort showed a time-dependent 
splenomegaly that started at week 3, while the group in- 
oculated with SiHac^y had no signs of splenomegaly at 
any time point (Figure 3C). Splenomegaly caused by 
SiHdiparentai tumors was associated with severe changes 
in the relative size of the red and white pulp and with 
infiltration of polymorphonuclear leukocytes in the 
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Figure 3 Growth of SiHacD\/ compared to S'lHaparentai a 
xenograft model in athiymic nude mice. (A) Tumor growth was 
determined in nu/nu mice bearing SiHacoi/ or SiHdparentai grown as 
xenograft at different time points post-inoculation of tine cells. Data 
were collected from two independent experiments (including a total 
of 15 mice/group) over a period of 5 weeks and are presented as 
the average tumor volume (mm^) ± SEM. Tumors were measured by 
means of a digital caliper in two directions (perpendicular diameters) 
and the formula V= (4nabV3, with 'a' and 'b' being the largest and 
smallest radius of the tumor, respectively, was applied to calculate 
the tumor volume. (B) Body weight of the mice (excluding the 
tumor mass) was determined after 5 weeks of inoculation of the 
cells and body weight gain was calculated as percentage of the 
initial body weight of the mice at the start of the experiment. Data 
shown represent the average (±SD) of five individual mice per 
group. Healthy control mice gained 47.0 + 4.4% of their initial body 
weight, while mice bearing the Si^dparentai cells only gained 33.2 ± 
7.1% (p < 0.01 compared to tumor-free mice). An intermediate body 
weight gain of 38.2 ± 9.8 (p > 0.05) compared to healthy control 
mice was recorded for mice with SiHacoi/- (Q Weight of the spleens 
was assessed at different weeks post-inoculation of the cells. 



extended red pulp (Figure 4). In the SiH^lcdv cohort, no 
considerable alterations in the morphology of the spleen 
were noted. 

When the different immune cell types in the spleens 
were quantified, striking differences were noted between 
animals inoculated with the two SiHa cell lines. Hence, 
the SiHsiparentai gfoup had a very pronounced and time- 
dependent increase in the number of neutrophils per 
spleen compared to healthy control mice (Figure 5). This 
increase could be attributed not only to splenomegaly 
but also to a higher percentage of neutrophils within the 
splenocytes in SiHsiparentai mice compared to healthy an- 
imals (data not shown). Splenomegaly in mice inoculated 
with SiWdiparentai was also associated with increased 
numbers of macrophages, NK-cells and B-cells per 
spleen compared to healthy animals. 

When comparing the SiHaci^y group with healthy 
animals, a slight increase (although significant) in the 
amounts of macrophages and neutrophils was measured, 
respectively, from weeks 3 and 4 onwards, while the 
numbers of NK-cells and B-cells did not differ between 
the two cohorts throughout the entire experiment 
(Figure 5). 

The SiHac^y cohort had a much lower amount of neu- 
trophils per spleen compared to the SiHsiparentai group 
which became evident from 2 weeks onwards. The num- 
ber of macrophages, NK cells and B cells were lower in 
animals harboring SiH^lcdv xenografts, starting from 
week 3 onwards, compared to the SiHdiparentai cohort. 

The groups with SiHsiparentai and SiHac^y tumors also 
differed in their ability to modulate tumor-derived 
cytokines and host-derived cytokines in the sera of 
mice (Figure 6). Most of the human-derived cytokines 
measured in the sera of tumor-bearing mice were 
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Table 3 Complete blood count of tumor-bearing mice 



Blood test 


Average ± SD 




p-value (t-test) 








Tumor- 
free 


SlHQparental 

tumor 


S'lHacDv 
tumor 


Tumor-free vs 

SlHaparental 


Tumor-free vs 
S'lHacDv 


S'lHSparental 

SIHacDv 


Hematology 














Hemoglobin [g/dl] 


1 3 g _l_ Q 3 




12.9±1.2 




n.s. 


n.s. 


Hematocrit 


U.^D X U.U 1 


U.j/ X U.UD 


U.^j X U.Uj 


n.s. 


n.s. 


n.s. 


HdL count LIU /ij 


p Qi _|_ n on 


7 OR + 1 17 
/ .UO X 1 . 1 / 


O.JO ± u.y 1 




n.s. 


n.s. 


rfn 


J 1 .Z X 1 .J 


DZ.y X U.J 


CI 7 -1- 1 Q 

J 1 ./ ± 1 .y 


n.s. 


n.s. 


n.s. 


IVILH Lpgj 


1 J.D ± U.J 


1 J 


1 j.y ± U.z 




n.s. 




IVILHL Lg/dlj 


30.6 ± 0.7 


28.2 ± 0.5 


29.8 ± 0.7 




n.s. 




Rn\A/ ro/ni 


1 7 4 _i_ Q g 


1 O.^ X U.O 


18 2 + 1 0 


n.s. 


n.s. 


n.s. 


Reticulocyte count 
[%] 


D. 1 X U.O 


Z.O X 1 .o 


J.J X U.J 


n.s. 


n.s. 


n.s. 


IRF [%] 


41.0 ±10.3 


24.8 ± 10.3 


36.7 ±5.4 


n.s. 


n.s. 


n.s. 


Platelet count [lOVl] 


1037± 113 


1566 ±22 


1 783 ± 240 






n.s. 


MPV [fl] 


7.1 ±0.2 


7.8 ± 0.4 


7.2 ±0.1 




n.s. 




WBC count [lOVl] 


2.1 7 ±0.22 


208.11 ±107.15 


9.39 ± 6.63 




n.s. 




Chemistry 














Creatinine [mg/dl] 


< 0.20 


< 0.20 


< 0.20 








AST [U/l] 


88±13 


392 ±188 


127±39 


•X 


n.s. 




ALT [U/l] 


48±11 


63 ±46 


33±16 


n.s. 


n.s. 


n.s. 


GGT [U/l] 


< 3 


12±6 


< 3 








Bilirubin [mg/dl] 


< 10 


< 10 


< 10 








CK [U/l] 


144± 16 


400 ± 60 


224 ±124 


**-x 


n.s. 


n.s. 


LDH [U/l] 


1263 ±361 


8966 ±1613 


3253 ± 1772 




n.s. 





Five weeks after inoculation of the tumor cells, blood was collected and examined by complete blood count. RBC: red blood cell; MCV: mean cell volume; MCH: 
mean corpuscular hemoglobin; MCHC: mean corpuscular hemoglobin concentration; RDW: red blood cell distribution width; IRF: immature reticulocyte count; 
MPV: mean platelet volume; WBC: white blood cell; AST: aspartate aminotransferase; ALT: alanine aminotransferase; GGT: gamma-glutamyl transferase; CK: creatine 
kinase; LDH: lactate dehydrogenase. Statistical significance was assessed based on unpaired two-tailed Student's t-test, n.s.: not significant; significant (p values 
< 0.05); very significant (p < 0.01); and ^^^: extremely significant (p < 0.001). 



undetectable. However, four tumor- derived cytokines 
(IL-6, IL-8, TNF-a, and IFN-y) were found time- 
dependently induced in the SiHsiparentai group. Human 
TNF-a and INF-y were undetectable at any time point 
in the SiHac^y group while human IL-6 and human IL-8 
were detected at very low levels starting at week 3. Not- 
ably, SiHsiparentai cells Were able to produce very high 
levels of IL-6. The decreased production of IL-6 by 
SiHdicDv compared to SiHdiparentai confirmed by an 
ELISA assay performed with cell-culture supernatants 
from both cell types (Additional file 6). 

Figure 6 also shows the mouse-derived cytokines 
(IL-lp, IFN-y, TNF-a, and IL-17) differently detected in 
tumor-bearing mice versus control healthy animals. In 
the tumor-free group, the host derived cytokines IL-1|3, 
IFN-y and TNF-a showed a time-dependent decline in 
concentration as a consequence of immune cell matur- 
ation while IL-17 levels remained stable in function of 
time. Sera of older mice showed lower and stable con- 
centrations of cytokines (data not shown), pointing to an 



effect of immune maturation on cytokine levels rather 
than an effect of housing or experimental conditions. 
Differences between healthy mice and the SiHsiparentai 
cohort started to be noticeable from week 2 (IL-lp), 4 
(TNF-a and IL-17) or 5 (IFN-y). While IFN-y, TNF-a 
and IL-17 levels remained stable or decreased during 
the first weeks in the SiHsiparentai group, their levels 
raised from week 4 onwards. Although IL-lp still de- 
creased in function of time in mice with SiHdiparentai 
tumors, the levels were significantly higher than in 
healthy animals. Except for IL-lp at week 5, host- 
derived cytokines levels did not significantly differ be- 
tween the tumor-free and the SiHac^y cohorts at any 
time point post-inoculation of the cells, pointing to a 
markedly diminished host inflammatory response com- 
pared to SiH2iparentai Xenografts. 

Discussion 

In the present study, we showed that SiHa cells that ac- 
quired CDV-resistance proved to be refractory to CDV- 
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Figure 4 Spleen pathology of mice inoculated with SiHacD\/ 
versus S'lHaparentai- Tissues were fixed in neutral buffered formalin, 
subsequently embedded in paraffin and 5 |jm sections were 
hematoxylin-eosin (H & E) stained and microscopically examined. 
Representative H&E stained tissue sections of spleens from tumor- 
free (A), and mice bearing SiHdparentai (B) or SiHacoi/ (Q xenografts, 
five weeks after inoculation of the cells (2.5x magnification). 
Splenomegaly caused by Sl^dparentai tumors was associated with 
severe changes in the relative size of the red and white pulp with 
an expanded red pulp, only sparse spots of white pulp and virtually 
no marginal zones remaining. In animals inoculated with SiHacoi// no 
major alterations in the morphology of the spleen were noted. 



antiproliferative effects and to CDV-induced apoptosis 
in vitro. These HPV-16 positive cervical carcinoma cells 
demonstrated a high barrier for the development of 
resistance to CDV as selection required prolonged 
exposure to CDV (approximately 45 passages during a 
2-years time). 

Genome wide gene expression analysis has been previ- 
ously used to identify gene expression signatures associ- 
ated with resistance to chemotherapeutic agents [30-32]. 
Here, we compared microarray gene expression values 
of SiHacDK with SiHsiparentai and bioinformatics analysis 
revealed the implication of a variety of biological func- 
tions and pathways (linked to cell death, cell growth and 
proliferation, cellular movement, metabolism, cell and 
tissue development as well as inflammatory response) 
changing following acquisition of resistance to CDV. 
Thus, it appears that acquisition of CDV-resistance is a 
multifactorial process, which is in agreement with find- 
ings on development of resistance to several chemother- 
apeutics [33]. 

By examining the identities of the genes in the 'in- 
flammatory response' exhibiting changes in expression 
upon acquisition of CDV-resistance, it can be assumed 
that the identified genes may not be the 'drivers' of 
drug-resistance, but they changed expression as a con- 
sequence of altered expression of the 'driver' genes. 
Candidate genes that should be further explored include 
c-Fos, c-Jun, PI3K and MAPK since they were changing 
expression upon acquisition of CDV-resistance and were 
involved in most of the inflammatory response pathways. 
Further investigations to elucidate the genes that drive 
acquisition of CDV-resistance are currently ongoing. 

The changes in inflammatory response observed in 
cells that acquire CDV-resistance are expected to be a 
consequence of the development of CDV-resistance ra- 
ther than the cause of the resistant phenotype in vitro. 
While not causing the resistant phenotype per se, the al- 
terations in inflammatory response are expected to affect 
the tumor microenvironment in vivo and to contribute 
to the observed reduction in pathogenicity and tumori- 
genicity. Therefore, it was interesting to investigate how 
acquisition of CDV-resistance in SiHa cells affected the 
inflammatory response induced by these cells in an athy- 
mic nude mice xenograft model. 

In the in vitro setting, SiH^lcdv proved clearly resistant 
to CDV but this must occur via a mechanism that does 
not directly involve cells of the immune system or the 
tumor microenvironment. In contrast, in vivo, the de- 
creased inflammatory response observed with SiHac^y 
compared to SiHsiparentai affected the tumor microenvir- 
onment and contributed to a reduced pathogenicity of 
the xenografts as SiHac^y provoked less inflammation 
in the xenograft model (evidenced by a reduced produc- 
tion of mice- and human-derived cytokines, diminished 
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Figure 5 Flow cytometric analysis of immune cell populations in the spleen of mice bearing SiHacD\/ or S'lHaparentai xenografts 
compared to healthy animals. Grl+CDl lb+ neutrophils, F4/80+CD1 lb+ macrophages, CD19+B220+ B-cells and CD3"DX5+ NK-cells are depicted 
as number of cells per spleen. Each bar represents the average (±SD) of three individual mice. 



effect on chemical and hematological blood parameters, 
lower number of immune cells in the spleen, and lesser 
splenomegaly compared to parental cells). 

In contrast to SiHac^vJ ^'^^^parentai generated a pro- 
nounced stimulation of immune cells (mostly neutro- 
phils) when evaluated in comparison to healthy animals. 
One could argue that the reduced induction of neutro- 
phils, macrophages, B-cells and NK-cells by SiHac^y 
could be the consequence of reduced growth rate ob- 
served for the SiHacDK not only in vitro but also in vivo. 
Yet, SiHdiparentai tumor size at week 3 was equivalent to 
that of SiHac^y at week 5 (tumor size of, respectively, 
351.6 ±259.8 mm^ and 342.0 ± 182.3 mm^) while the 
amount of neutrophils, macrophages and NK-cells was 
considerable higher in mice with SiHsiparentai xenografts 
than in those with SiH^lcdv tumors at these time points. 
Similarly, when putting side by side the SiHdiparentai 
SiHdiCDv groups at the moment that they have an 
equivalent tumor size (week 3 and week 5, respectively), 
IL-1|3 was detected in higher amounts in the SiHdiparentai 
cohort. IL-1|3 plays a key role in the regulation of 
neutrophil recruitment through up-regulation of endo- 
thelial adhesion molecule expression on endothelium 
and through induction of local chemokine production 
(including IL-8) production [34], and indeed lower IL-1|3 



levels correlated with lower numbers of neutrophils in 
the SiHac^y cohort. 

Neutrophils and macrophages have a major role in 
defense mechanisms and protect the host from injury 
and infections. However, they were shown to infiltrate 
most solid cancers and tumor-associated macrophages 
(TAMs) and tumor-associated neutrophils (TANs) were 
shown to be involved in stimulation of tumor growth, 
their densities being linked to poor outcomes and 
shorter survival in several cancer types [35,36]. A recent 
study showed that elevated white blood cells and neutro- 
phil counts at the time of recurrence diagnosis corre- 
lated with shorter survival in patients with recurrent 
cervical cancer [37]. In other cancers, such as colon can- 
cer, small cell lung carcinoma, and melanoma, an ele- 
vated neutrophil-to-lymphocyte ratio also predicted a 
significantly higher risk of death [38-40]. 

Recently, a role for the spleen as a site for storage and 
rapid deployment of monocytes to inflammatory sites 
has been unraveled, identifying splenic monocytes as a 
resource that the body uses to regulate inflammation 
[41]. Cortez-Retamoza and colleagues [42] demonstrated 
the function of the spleen as a reservoir of monocytes 
using a mouse model of lung adenocarcinoma. High 
numbers of TAMs and TANs relocated from the spleen 
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Figure 6 Host and human cytokine levels in the sera of mice harboring SiHacD\/ or S'lHaparentai tumors compared to tumor-free animals. 

Cytokine levels in the sera of mice were determined after 1-5 weeks post-inoculation of the cells. Each bar represents the average (+SD) of three 
individual mice. Human cytokines GM-CSF, IL-lb, IL-2, IL-4, IL-5, IL-7, IL-IO, IL-1 2, IL-1 3, IL-1 7, MIPip, M-CSF, TNF-p were undetectable. The mouse 
cytokines IL-2, IL-4, IL-5, IL-10, granulocyte-macrophage colony-stimulating factor (GM-CSF), IL-6, IL-1 3, macrophage colony-stimulating factor 
(M-CSF) and IL-1 8 were detectable but no differences among the three groups of mice were observed. 



to the tumor stroma. Furthermore, removal of the spleen 
(either before or after tumor initiation) reduced TAMs 
and TANs responses markedly and delayed tumor 
growth [42]. Local accumulation of granulocytes and 
macrophage progenitors in the splenic red pulp was 
linked to the reservoir capacity of the spleen during 
tumor progression. Our data showing an infiltration of 
polymorphonuclear leukocytes in the extended red pulp 
in the SiHsiparentai Xenograft cohort (but not in the 
SiHac^y one) suggest that the spleen might also play an 
important role as reservoir of monocytes. Moreover, a 
pronounced increase in the number of WBC was de- 
tected in the SiHsiparentai but not in the SiH^lcdv group. 
Our microarray data also indicated that acquisition 
of CDV-resistance was associated with reduction of 
'inflammatory response', 'activation of granulocytes', 
inflammation of organ' and activation of neutrophils' 
(Additional file 2), which can explain the diminished 
stimulation of the production of neutrophils and 
macrophages by the host. Decreased expression of 
genes whose products are responsible for activation 
of neutrophils and/or granulocytes (such as complement 
components, endothelin 1, IL-1 5, integrin p2, monocyte 
chemoattractant protein 1, macrophage inflammatory pro- 
tein 4a, protein kinase C inhibitor 2, G TP-binding protein 
Ram, and monocyte differentiation antigen CD 14) point 
to a decreased capacity of SiHacoy cells to activate and 



attract neutrophils and macrophages at the tumor site 

compared to SiH2iparental' 

Overafl, our data showed that SiH^lcdv elicited a re- 
duced inflammatory response in the xenograft model 
when evaluated in comparison with SiH^Lparentai- Inflam- 
mation is present in almost aU cancer tissues and the 
inflammatory state is necessary in tumor tissue remodel- 
ing, angiogenesis and metastasis [43-45]. Altered expres- 
sion of cytokines and growth factors is crucial in the 
malignant transformation of many cancers. Inflamma- 
tion, actually smoldering' inflammation, is now consid- 
ered as one of the hallmarks of cancer [43,46]. Recent 
studies pointed out the importance of cytokine profiles 
in patients with cervical intraepithelial and invasive neo- 
plasia, suggesting that tumor progression is dependent 
on suppression of cellular immunity [47,48]. Hence, de- 
creased levels of Thl cytokines were reported in high- 
grade lesions, consistent with the role of Thl cytokines as 
potent activators of cell-mediated immunity [48-50]. Scott 
and colleagues also demonstrated that persistence of an 
HPV infection is linked to a failure to express Thl cyto- 
kines [51]. Chronic Th2 type inflammation is commonly 
seen during persistent infection with high-risk HPV types 
promoting tumor progression [52] . Furthermore, high-risk 
HPV types are able to initiate a local Th2 inflammation at 
an early stage, creating an immunosuppressive microenvir- 
onment that contributes to tumor progression [47] . 
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We have previously shown that the production of 
a number of cytokines by SiHsiparentah including the 
pro-inflammatory cytokines IL-6, IL-8, TNF-a and IFN- 
y, is decreased following CDV therapy in the xenograft 
model in nu/nu mice [24]. Here, we demonstrated that 
SiHacDK produced significant lower levels of these pro- 
inflammatory cytokines in mice. These findings were 
supported by bioinformatics analysis of microarray gene 
expression profiling that showed alteration of interleukin 
(IL-1, IL-6, IL-8, IL-9, IL-10) and interferon signaling 
pathways. 

Acquisition of CDV-resistance resulted in inhibition of 
the IL-6, IL-9, and IL-10 signaling pathways as inferred 
by a decreased expression of STATS, SOCS2 and 
SOCS3. The STATS protein is activated through phos- 
phorylation in response to various cytokines and growth 
factors including IFNs, EGF, IL-5, and IL-6, mediating 
the expression of a variety of genes in response to cell 
stimuli, and thus playing a key role in many cellular 
processes [5S,54]. SOCS family members are cytokine- 
inducible negative regulators of cytokine receptor signal- 
ing via the Janus kinase/signal transducer and activation 
of transcription pathway (the JAK/STAT pathway) [55]. 
Transcripts encoding SOCS are upregulated in response 
to cytokine stimulation, and the corresponding SOCS 
proteins inhibit cytokine-induced signaling pathways. 
Therefore, SOCS proteins form part of a classical nega- 
tive feedback circuit [56,57]. Expression of SOCS2 can be 
induced by a subset of cytokines such as GM-CSF, IL-10 
and IFN-y while that of SOCSS by IL-6, IL-10 and IFN-y. 
It can be inferred that reduced expression of STATS and 
SOCS genes in SiHac^y versus SiHsiparentai is the conse- 
quence of reduced levels of cytokines, and indeed, 
SiHac^y produced lower levels of pro-inflammatory cyto- 
kines (IL-6, IL-8, TNF-a and IFN-y) in mice. 

In the xenograft model, human IL-6, IL-8, and TNF-a 
are expected to have an important role in the mice path- 
ology because they are known to be biologically active in 
mice, in contrast to IFN-y and its receptor that are spe- 
cies specific [58]. SiHsiparentah but not CDV-resistant 
cells, produced high levels of IL-6. This cytokine is 
known to induce extensive extramedullar hematopoiesis 
leading to production of neutrophils that localize to the 
tumor microenvironment promoting tumor growth by 
protease-induced angiogenesis [59]. 

TNF, originally identified for its ability to induce rapid 
hemorrhagic necrosis of experimental tumors, is now 
recognized as a central mediator of inflammation, 
representing one of the molecular links between chronic 
inflammation and the subsequent development of malig- 
nant disease [60]. TNF-a is a strong activator of NF-kB, 
an injury transcription factor that contributes to cell sur- 
vival, proliferation, invasion, inflammation and angio- 
genesis [61]. Tumor promotion by TNF-a can involve 



diverse pathways, including enhancement of tumor 
growth and invasion, leukocyte recruitment, angiogen- 
esis and facilitation of mesenchymal transition [62]. 
SiHac^y showed increased expression of the TNF recep- 
tor TNFRSFllB and diminished expression of the TNF 
ligand TNFSF15 (Additional file 4), which is expected to 
affect NF-kB activation and apoptosis induction. This 
hypothesis is based on the fact that TNFRSFllB is a 
decoy receptor for RANKL (receptor activator of NF-kB 
ligand) and TRAIL (TNF-related apoptosis-inducing lig- 
and), and that TNFSF15 (which is inducible by TNF and 
IL-la) binds to TNFRSF21 (an activator of NF-kB and 
of apoptosis). Further evidence for an effect on NF-kB 
activation in SiHac^y versus SiHsiparentai is provided by 
increased expression of the TNF associated factor 
TRAFS (a known inhibitor of NF-kB activation) and of 
IKBKG [the regulatory subunit of the inhibitor of kappa 
B kinase (IKK) complex, also known as NEMO]. 

The decreased expression of several genes implicated 
in the HMGBl (high mobility group box 1) signaling 
pathways in SiH^cDv versus SiHsiparentai further supports 
the reduced tumorigenicity and inflammation of cells that 
acquired CDV-resistance. As post-translational modifi- 
cations determine intracellular distribution and key 
functions of HMGBl, changes at the mRNA level for 
HMGBl were not detected. However, in the HMGBl 
signaling pathway, expression of mitogen-activated pro- 
tein kinases (MAPKs) and of the serine/threonine kin- 
ase AKTS was reduced in SiH^lcdv versus SiHsiparentah 
leading, respectively, to diminished expression of c-Fos 
and c-Jun and to regulation of NK-kB. c-Fos and c-Jun 
form the transcription factor complex AP-1 which regu- 
lates gene expression in response to a variety of stimuli 
(such as cytokines, growth factors, stress, and microbial 
infections) and controls a number of cellular processes. 
HMGBl, considered as a prototypic damage-associated 
molecular pattern (DAMP) molecule, acts as both a 
ligand and a sensor of the signal-transducing innate 
responses. Therefore, it can be assumed that a decrease in 
HMGBl signaling following acquisition of CDV-resistance 
may result in lower stimulation of pro-inflammatory 
cytokines. 

Another interesting finding when comparing SiHac^y 
and SiHdiparentai is their differences in TLR signaling, 
with TLRS and TLR4 is downregulated in SiH^lcdV' 
TLRs activate several signaling elements that results in 
activation of pro -inflammatory cytokines, regulating 
apoptosis, antimicrobial response and immune respon- 
ses. Expression of TLRs in tumor cells can promote 
inflammation and cell survival in the tumor micro- 
environment [6S,64]. Moreover, expression of TLRs in 
esophageal squamous carcinoma [65] and in cervical le- 
sions [66] was shown to correlate with disease severity. 
As TLRs promote tumor cell growth and cytokine 
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secretion, leading to the escape of tumor cells from im- 
mune surveillance, it can be assumed that reduced TLR 
expression in SiH^lcdv will contribute to a reduced in- 
flammatory response and decreased tumor growth com- 
pared to the parental cells. 

Further evidence for lower tumorigenicity induced by 
SiHac^y versus SiHsiparentai in mice is provided by 
changes in the 'MSP/RON signaling pathway . Macro- 
phage-stimulating protein (MSP) activates the RON re- 
ceptor tyrosine kinase, which regulates several activities 
of epithelial cells [67]. The MSP-RON pathway plays 
also a role in epithelial carcinogenesis and RON is found 
over-expressed in many breast, colon, and pancreatic 
tumors [67]. As activation of the MSP-RON pathway 
directs invasive growth (characterized by increased cell 
replication, migration, and matrix invasion) [68,69], it 
can be inferred that the decreased expression of genes 
involved in this pathway [such as TLR4, TLR3, monocyte 
chemoattractant protein 1 (CCL2) and integrin |32] in 
SiHsiCDv versus SiHsiparentai will be translated in a re- 
duced tumorigenicity in vivo. 

In the context of developing CDV as an anti-cancer 
drug, our findings have therapeutic/biological signifi- 
cance since we showed that acquisition of CDV- 
resistance is expected to result in a reduced malignant 
phenotype. 

Today, no evidence for the development of resistance 
to CDV in the treatment of HPV-associated (malignant) 
lesions has been reported. 

Conclusions 

Although several studies have characterized CDV-resistant 
herpes- [70], and poxviruses [71], this is the first study 
reporting the in vivo characterization of tumor cells 
selected for CDV-resistance. Similarly to a reduced 
pathogenicity described for CDV-resistant viruses, de- 
velopment of resistance to CDV as an anti-cancer 
drug was associated with a marked reduction in 
pathogenicity. The present study contributes to our 
understanding on how the alterations in inflammatory 
response following acquisition of CDV-resistance while 
not causing the resistant phenotype per se affect the tumor 
microenvironment in vivo and contribute to a reduced 
pathogenicity and tumorigenicity. 

Methods 

Compounds 

Cidofovir [CDV, (5)-HPMPC, (.S)-l-(3-hydroxy-2-phos- 
phonylmethoxypropyl)cytosine], PMEG {[9-(2-phospho 
nylmethoxyethyl)guanine]} and cPr-PMEDAP {cyclo- 
propyl-9- [2-(phosphonomethoxyethyl] diaminopurine] } 
were kindly provided by Gilead Sciences, Inc., Foster City, 
California. Cytarabine [Ara-C, (|3-D-Arabinofuranosyl) 
cytosine] was obtained from Sigma. 



Cells 

SiHa cells, HPV16-positive cervical carcinoma (ATCC, 
# HTB-35™), were maintained in Dulbeccos modified 
Eagle s medium supplemented with 10% fetal calf serum. 
SiHa cells resistant to CDV were selected by passing the 
cells under increasing drug concentration for approxi- 
mately 45 passages (initial drug concentration of 1.6 [iM 
and final concentration of 317.3 [iM) during a 2-years 
time. The parental SiHa cells and those selected for re- 
sistance to CDV were denoted SiHsiparentai and SiHac^i^ 
respectively. In order to demonstrate that both cell lines 
were related, short tandem repeat (STR) analysis was 
performed at the Forensic Laboratory of UZ Leuven 
(Leuven, Belgium). Despite some small alterations fol- 
lowing long-term culturing of the cells, the STR analysis 
confirmed that SiHsiparentai and SiH^lcdv were related 
and thus that the resistant cell line is indeed a derivative 
of the parental cell line (see Additional file 7). 

Drug-antlproliferative effects and in vitro growth rate 

Inhibition of SiHsiparentai and SiH^lcdv growth was deter- 
mined following different times of incubation with the 
compounds. Compounds were tested at different con- 
centrations in a range of 0.63 [xM - 634.7 [xM for CDV, 
0.0065 [iM - 6.54 [iM for PMEG, 0.061 [iM - 60.99 [iM 
for cPr-PMEDAP, and 0.0205 [iM - 20.53 [iM for Ara-C. 
Antiproliferative effects were expressed as CC50 (50% 
cystostatic concentration), or concentration required to 
reduce cell growth by 50%. 

Doubling time (DT) of SiHsiparentai and SiH^lcdv was 
determined in 48-well microtiter plates from growth 
curves performed in absence of the drug by using the 
formula: DT = (t - to)/(log2N - log2No), where t and to 
are the times at which the cells were counted, and N 
and No are the cell numbers at times t and to. 

Detection of apoptosis 

To differentiate between living, apoptotic and necrotic 
cells, SiHdiparentai and SiHdicDv were grown for 7 days in 
the presence of CDV or PMEG. Cells were simultan- 
eously stained with annexin V-FITC and propidium iod- 
ide (PI) using the FITC Annexin V Apoptosis Detection 
Kit (BD Pharmigen™). Dual fluorescence dot plots 
were determined with a FACSCalibur flow cytometer 
equipped with CellQuest software (BD Biosciences). 

Microarray experiments 

SiHdiparentai and SiHsicDv cells were allowed to grow for 
72 h in medium without CDV. Total RNA of 1 x 10^ 
cells was isolated with TRIzol reagent (Invitrogen) ac- 
cording to the manufacturers instructions. The RNA 
was further purified by RNeasy Mini Kit (Qiagen). RNA 
quality and quantity were assessed by using a Bioanaly- 
zer system (Agilent). 
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Human Genome U133 Plus 2.0 arrays (Affymetrix) 
containing more than 54,000 probe sets and covering 
approximately 38,500 genes were used to analyze the ex- 
pression profile of the two cell lines, and both conditions 
were tested in triplicate. Array hybridization, scanning 
and image analyzing were done following the manufac- 
turer s protocols (Affymetrix GeneChip Expression Assay) 
at the VIB Nucleomics Core Facility (www.nucleomics.be). 
Raw data were corrected for background signal using the 
RMA (Robust Multi- array analysis) algorithm (affy_1.22.0 
package of BioConductor). The detection (Present/ 
Absent) call generated by the Affymetrix microarray 
suite version 5 software (MAS 5.0) was used to re- 
move probe sets that were not reliable detected in any 
of the microarrays before further analysis. 

Differentially expressed (DE) probe sets between 
SiWdiparentai ^nd SiHac^y were determined using a mod- 
erated t-statistic test [LIMMA (linear models for micro- 
array data), BioConductor]. The Benjamini-Hochberg 
correction for multiple testing [p < 0.05, false discovery 
rate (FDR) = 0.05] was performed. Probe sets were 
considered significantly DE if the absolute fold-change 
(FC) was > 2 and the P-value was < 0.05 (LIMMA) after 
applying the Benjamini-Hochberg correction. 

The entire set of microarray data is deposited in the 
Gene Expression Omnibus (GEO, http://www.ncbi.nlm. 
nih.gov/projects/geo) according to MIAME standards 
under accession number GSE26748: http://www.ncbi. 
nlm.nih.gov/geo/query/ acc.cgi?token=lpivfquymowyazo& 
acc=GSE26748. 

Bioinformatics analysis of differentially expressed genes 
was carried out with Ingenuity Pathways Analysis (IPA, 
Ingenuity® Systems) version 9. Data sets with the corre- 
sponding FC and P-value were uploaded into the IPA 
(Ingenuity Pathway Analysis, Ingenuity® Systems) soft- 
ware. Stringent criteria, equivalent to those described 
for the selection of DE probes, were applied to identify 
DE genes. When genes were represented by 2 or more 
probe sets on the arrays, only the maximum FC was 
used. Uncharacterized probe sets were not included in 
the analysis. 

The IPA application reveals relevant pathways and bio- 
logical functions by comparing the number of genes that 
participate in a given function or pathway, relative to the 
total number of occurrences of those genes in all the 
pathways stored in the IPKB (Ingenuity Pathway Know- 
ledge Base). Validation of the microarray data was per- 
formed with 4 genes (randomly selected) by quantitative 
reverse transcription-polymerase chain reaction (qPCR) 
as previously reported [28]. 

Animal experiments 

Female nu/nu NMRI mice (4-5 weeks old) were pur- 
chased from Janvier Breeding Center. All animal work 



was approved by the KU Leuven Ethics Committee for 
Animal Care and Use (Permission number: P160-2008). 

Mice were inoculated sub-cutaneously on the back 
with 2 X 10^ cells in a volume of 200 (il, week 0 being 
considered the time point of cell inoculation. To esti- 
mate body weight gain, mice were sacrificed weekly and 
tumors were excised, weighed and subtracted from the 
total body weight. Gain in body weight was calculated as 
the percentage of body weight gained compared to the 
mice weight at week 0. 

Spleens from 3 mice per group were isolated at differ- 
ent weeks to determine the percentage of immune cell 
populations. Spleens were processed and splenocytes 
were stained with specific antibodies and analyzed by 
flow cytometry as described previously [24]. One mouse 
per group was euthanized weekly to collect various tis- 
sues for histopathological examination. 

Total blood from 5 mice per group was collected in 
EDTA tubes at week 5 to perform hematology and blood 
chemistry testing at the University Hospitals Leuven, 
Department of Laboratory Medicine, Leuven, Belgium. 

At various time points, tumor- and host-derived cyto- 
kines were quantified in the sera of mice (3 animals per 
group) with a Bio-plex 200 system (Bio-Rad Laboratories) 
according to the manufacturer s protocols. 

Statistical analysis 

Statistical significance was assessed based on unpaired 
two-tailed Students ^-test with GraphPad Prism 5 soft- 
ware (GraphPad Software Inc., La Jolla, CA, USA). Sig- 
nificance was indicated as: ns, not significant (p > 0.05); 

significant (p < values 0.05); very significant (p < 
0.01); and extremely significant (p < 0.001). 

Additional files 



Additional file 1: Validation of gene expression between 
microarray and qPCR. The gene expression levels of microarray are 
presented by log2 fold changes, whereas those of qPCR are indicated by 
AACt that are comparable to the log2 fold change values in microarray. 
DHRS2 (Dehydrogenase/reductase SDR family member 2), HIST1H2AC 
(Histone H2A type 1-C), GALC (Galactocerebrosidase), KRT14 (Keratin 14). 

Additional file 2: Functional annotations associated with 
inflammatory response found to be distinct between S\Wdiparentai 
and SiHacD\/. The criteria for selection of functional annotations were 
based on z-score and statistical significance (P-value < 0.05). The regula- 
tion z-score predicts whether an identified biological function is activated 
or inhibited. Positive z-scores indicate activation of a biological function, 
while negative z-scores suggest an inhibition. Absolute z-score values 
above 1 were considered significant. 

Additional file 3: Canonical pathways related to immune response 
found to be distinct between SiHdiparentai and SiHacoi/. The 

significance of the associations between the genes from the two data 
sets and the canonical pathways were determined based on two 
parameters: (a) the P-value, calculated by the Fischer's exact test, that 
determines the probability that there is an association between the 
genes in the data set and the canonical pathway that cannot be 
explained by chance alone and (b) the ratio of the number of genes 
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from the data set in a given pathway divided by the total number of 
molecules in the given canonical pathway. P-values < 0.05 were 
considered statistically significant. 

Additional file 4: Gene expression changes related to the 
'inflammatory response' function in SiHacD\/ compared to 
S'lHa parental' Genes were considered significantly differentially expressed if 
the absolute fold-change (FC) was > 2 and the P-value was < 0.05 
(LIMMA) after applying the Benjamini-Hochberg correction. Upregulated 
and downregulated genes are indicated by respectively positive and 
negative log2 fold changes. 

Additional file 5: Effects of intratumoral CDV treatment on tumor 
growth in xenograft model. (A) Kinetics of tumor growth were 
determined in mice bearing a xenograft that received either no 
treatment or intratumoral CDV treatment (25 mI of a 10 mg/ml CDV 
solution) once a day, five times per week, for a period of four weeks. 
Data are presented as the average tumor volume [mm^] of five mice 
(±SEM). Tumors were measured by means of a digital caliper in two 
directions (perpendicular diameters) and the formula V = (4nab'^)/3, with 
'a' and 'b' being the largest and smallest radius of the tumor, respectively, 
was applied to calculate the tumor volume. (B) Efficacy of treatment 
following four weeks of intratumoral CDV treatment. Efficacy was 
evaluated by means of the inhibitory rate (IR) on tumor growth and was 
calculated as IR = (C-T)/C x 100% (C and T being the respective tumor 
volumes of the untreated controls and the treated tumors). Treatment of 
SiUdparentai xenografts resulted in an IR of 95% (p-value < 0.01), while 
treatment of Si Hacoi/ xenograft showed an IR of 51% (p-value > 0.05). 

Additional file 6: Interleukin 6 (IL-6) levels in S'lHaparentai and 
SiHacD\/cell culture supernatants. IL-6 levels were measured using the 
enzyme-linked immunosorbent assay kit (Invitrogen™) using cell culture 
supernatants. Cells were seeded at a density of 4 x 10"^ cells per well in 
3 ml culture medium in 6-well plates. After 24 h the medium was chan- 
ged, and supernatant was subsequently collected after 3, 5, and 7 days. 
Human IL-6 levels in the supernatant (diluted 1:25 in PBS) of cultured 
cells were determined by using the IL-6 human ELISA kit following 
manufacturer's instructions. Samples were measured in triplicate. 

Additional file 7: STR profile of SiHapo^-enfo/ and SiHacoi/ cells. The 

different alleles for the STR loci that were identified in both cells lines are 
shown in the table. Determination of the STR profile of both cell lines 
illustrated a driftof a few markers (i.e. D1S1656, D21S11 and D1S1677) 
following long-term culturing of the cells. Overall these data demon- 
strated the relationship between the two cell lines and confirmed 
SiUdcDv being a derivative of the SiUaparemai cells. 
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